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»^, We study the Sivers effect in transverse single spin asymmetries (SSA) for pion and kaon 

C^ ; production in Semi-Inclusive Deep Inelastic Scattering (SIDIS) processes. We perform 

(D ; a fit of A'jjy''^''^'^^' taking into account the recent data from HERMES and COMPASS 

1-^ ■ Collaborations, which allow a new determination of the Sivers distribution functions for 

quark and anti-quark with u, d and also s flavours. Estimates for forthcoming SIDIS 

experiments at COMPASS and JLab are given. 

>, 

^^ ' Data on the transverse single spin asymmetry A^'j! for polarized SIDIS processes, 

^^ ■ £p{S) ^ £'hX, collected by the HERMES [J and COMPASS [2] Collaborations allowed 

us [3111] to perform a rather weU constrained extraction of the Sivers distribution function 
[SJ [5] for u and d quarks, assuming a negligibly small Sivers sea. Recently, much higher 
f^ ' statistics data on the A^'j^ ' azimuthal asymmetries for SIDIS have become available: 

OO . in Ref. [7] the HERMES Collaboration presents neutral pion and charged kaon azimuthal 

asymmetries, in addition to higher precision data on charged pion asymmetries; moreover, 
Refs. [SI [S] show the COMPASS Collaboration measurements for separated charged pion 
^.j^ ■ and kaon asymmetries, together with some data for Kg production. 

^ ■ Here we present the analysis of these new experimental sets of data [TO] ■ They give us a 

P^ \ better understanding of the u and d flavour Sivers distribution functions at low-intermediate 

X and, most importantly, a first insight into the sea contributions to the Sivers functions. 

The SIDIS transverse single spin asymmetry A^^'j^ ' is defined as 

.sin(0^-0s) _r, ! d(j)s d(j)h [da^ - da^] sin(0,t - (j)s) 

where (ps and 4>h are the azimuthal angles identifying the directions of the proton spin S 
and of the outgoing hadron h in the 7*p cm. frame, see Fig. 1 of Ref. [IQ^. Taking into 
account intrinsic parton motion, this transverse single spin asymmetry, can be written, at 
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(2) 
where Lp defines the direction of the incoming (and outgoing) quark transverse momentum, 
fc^ = /ex(cosiy9, sin(y9, 0); fq/p{x,kj_) is the unpolarized x and k± dependent parton distri- 
bution function (PDF); da^'^^^'^ / dQ^ is the unpolarized cross section for the elementary 
scattering £q -^ £q; D^{z^pi_) is the fragmentation function describing the hadronization of 
the final quark q into the detected hadron h with a light-cone momentum fraction z and a 
transverse momentum p^ with respect to the fragmenting quark; finally, A^/^/pt {x, k±) is 
the Sivers function, parameterized in terms of the unpolarized distribution function as: 

A%p, {x, k^) = 2Nq{x) h{k^) fq/p{x, k^) , (3) 



with 



Nqix) = N.x-^il - xf^ ("-+f:^p"'^^ , Kk,) ^ ^/2^^ e-i/-? , (4) 



iq h'q 



where Nq £ [—1, 1], aq, f3q and Mi (GeV/c) are free parameters to be determined by fitting 
the experimental data. Notice that h(k±) < 1 for any k± and |A/'g(.'i;)| < 1 for any x, therefore 
the positivity bound for the Sivers function is automatically fulfilled. For the unpolarized 
distribution and fragmentation functions, we adopt the common factorized gaussian form 

fq/p{x, k^) = fq{x) -i^ e-'=i/('=l) , D';,{z,p^) = D^q{z) -^^ c-^/'^pI) , (5) 

with {k\) = 0.25 (GeV/c)2 and {p\) = 0.20 (GeV/c)^ fixed by analysing the Cahn effect 
in unpolarized SIDIS, as in Ref. [Sj. The parton distribution functions fq{x) and the frag- 
mentation functions Dq{z) also depend on Q^ via the usual QCD evolution, which will be 
taken into account, at LO, in all our computations. 

Fragmentation functions are a crucial ingredient of our fit. We have considered three 
different sets: KRE [n], HKNS [12] and DSS [13]. All these sets are basically equivalent 
as far as pion asymmetries are concerned. However there are important differences in the 
description of kaon data. In particular the DSS set, contrary to the other two sets, is such 
that D^ {z) ^ D!^ (z) over the whole z range. This feature is crucial when studying kaon 
production processes: first of all, it allows to reproduce kaon multiplicities at HERMES; 
secondly, it enables us to achieve kaon asymmetries larger than those corresponding to pion 
production. For these reasons we have chosen the DSS set for our fit. Contrary to the 
fragmentation sector, the use of different sets of unpolarized distribution functions does not 
affect our results significantly; here we use the GRV98LO set [14] . As the SIDIS data 
from HERMES and COMPASS have a limited coverage in x, typically a; < 0.3 - 0.4, the 
experimental asymmetries we are fitting contain very little information on the large x tail 
of the Sivers functions. Therefore we assume the same value of /3 (which is related to the 
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Table 1: Best values of the free parameters for the 'broken sea' ansatz. The errors are 
determined according to the procedure explained in Appendix A of Ref . [10] . 
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shape of the Sivers functions at large a;) for all Sivcrs functions, setting psea ~ f3u ~ Pd = P- 
Notice that this choice artificially reduces the width of the uncertainty band at large x. 
Moreover we assume the same a = a^ea for all sea quarks. Thus for this so called 'broken 
sea' ansatz fit we then have 11 parameters. The results we obtain for these parameters 
by fitting simultaneously the four experimental data sets on A^'^ , corresponding to 

pion and kaon production at HERMES [7] and COMPASS ID, are presented in Table I 
together with the correspond- 
ing errors, estimated accord- 
ing to the procedure outlined 
in Appendix A of Ref. [TU] . 

The fit performed under 
the 'broken sea' ansatz shows 
a good description of pion and 
kaon asymmetries. We ob- 
tained x^ = 1-20 per data itj . itj 
point for A'+ production at 
HERMES [7], while for pi- 
ous we had x^ = 0-94 per 
data point, for a total Xdof ~ 
1.00. Our results confirm that 
A^/„/pT > and A^fd/pr < 
as found in Ref. [J. More- 
over HERMES data on kaon 
asymmetries cannot be ex- 
plained without a sea-quark 
Sivcrs distribution. In partic- 
ular we find that A^/j/pT > 
0. Using the Sivers func- 
tions determined through our 
fit, we have given predictions 
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Figure 1: The Sivers distribution functions for u, d and 
s flavours as determined by our simultaneous fit of HER- 
MES and COMPASS data. On the left panel, the first mo- 
ment a; A^/(i)(a;) = ~x f^r^'^^'' {x){ see Eq. 17 of Ref. [TO]) 
is shown as a function of x for each flavour, as indi- 
cated. Similarly, on the right panel, the Sivers distribution 
x A'^f{x,k±) is shown as a function of k± at a fixed value 
of X for each flavour, as indicated. 



for COMPASS 



experiment operating with a 

hydrogen target and at JLab, on proton, neutron and deuteron transversely polarized tar- 
gets; for details see Ref. [TO] . 

We have performed a comprehensive analysis of SIDIS data on Sivcrs azimuthal de- 
pendences. It turns out that the data, and in particular the unexpectedly large value of 
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^sirH<ph rps) £qj, j^+^ demand a non vanishing, and large, Sivers distribution for s quarks. 
The other sea quark (u, d, s) contributions are less well determined, although they also seem 
to be non vanishing. 
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